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ABSTRACT 

We present high-resolution CARMA 230 GHz continuum imaging of nine deeply embedded proto- 
stars in the Serpens Molecular Cloud, including six of the nine known Class protostars in Serpens. 
This work is part of a program to characterize disk and envelope properties for a complete sample 
of Class protostars in nearby low-mass star forming regions. Here we present CARMA maps and 
visibility amplitudes as a function of ^-distance for the Serpens sample. Observations are made in 
the B, C, D, and E antenna configurations, with B configuration observations utilizing the CARMA 
Paired Antenna Calibration System. Combining data from multiple configurations provides excellent 
^-coverage (4 — 500 kA), allowing us to trace spatial scales from 10 2 to 10 4 AU. We find evidence 
for compact disk components in all of the observed Class protostars, suggesting that disks form at 
very early times (t < 0.2 Myr) in Serpens. We make a first estimate of disk masses using the flux at 
50 kA, where the contribution from the envelope should be negligible, assuming an unresolved disk. 
The resulting disk masses range from 0.04 M to 1.7 M , with a mean of approximately 0.2 Mq. Our 
high resolution maps are also sensitive to binary or multiple sources with separations > 250 AU, but 
significant evidence of multiplicity on scales < 2000 AU is seen in only one source. 
Subject headings: stars: formation — ISM: individual (Serpens) — submillimeter — techniques: in- 
terferometric 



1. INTRODUCTION 

Circum-stellar and circum-protostellar disks are an im- 
portant component of the star formation process, trans- 
ferring mass from the surrounding envelope to the grow- 
ing central protostar. While disk formation is a natu- 
ral result of collapse in a rotating core, it is still not 
known how soon after protostellar formation the disk 
appears, or how massive it is at early times. Simple an- 
alytical theories of star formation suggest that centrifu- 
gally supported disks should start out quite small (radius 
< 10 AU), and consequently with very low mas, and grow 
with time (jTerebev. Shu, fc Cassenlll984f ). More recent 
simulations of self-gravitating and viscous disks indicate 
that relatively massi ve disks can form at very early times 
(e.g. I Vorobvovl 120091 ). Magnetically support ed disks can 
also be much larger (radii up to 1000 AU; iGalli fc Shul 
1993(), and thus more massive at early times. 

To understand the formation and growth of disks, and 
their relationship to the accretion process, it is critical 
to observe and characterize disks in embedded sources, 
where most of the mass of the growing star is being ac- 
creted. Directly observing disks in the embedded phases 
is much more difficult than in more evolved pre-main 
sequence stars, as they are hidden within the dense, ex- 



tinct ing protostellar envelopes. 

Embedded protostars are defined here as protostars 
with substantial (> 0.1 Mq) envelope masses remain- 
ing. We will also refer to young so urces by class, where 
the bolometric te mperature T bo ^ (jMvers fc La dd 1993; 
iChen et aLl [T995) is used to divide sources into Class 
(Tboi < 70 K), Class I (70 K< T boi < 650 K), and 
Class II (650 K< T boi < 2800 K). Here we assume that 
these classes correspond to an evolu tionary sequence (e.g. 
iRobitaille etaLl 120061 : 1 Andrei 1 19941 ): in Class the pro- 
tostar has accreted less than half its final mass (M* < 
M env \ in Class I more than half of the mass has been 
accreted (M* > M env \ and in Class II the envelope has 
dispersed, leaving only a circum-stellar disk. 

In the last two decades a number of studies 
have focused on detecting disks in embedded pro- 
tostars (e.g. [Chandler et all 119951: iBrown et al.l 1 2000: 
Loonev et al. 2000; Harve vet al.l 1 2003: J^rgen sen et al.l 
2005bl : lEisner et all 120051 : lAndrews fc Williamsl 120071 ). 

1 The bolometric temperature is the temperature of a black- 
body with the same mean frequency as the observe d spectral e nergy 
distribution (T boi = 1.25 x lCT 11 ^) K: Myers fe Ladd 1993) and 
is a better me asure of the evolut ionary state than the infrared 
spectral index ([Enoch et al.H2009al V 
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TABLE 1 

Serpens Targets 



ID 


c2d name/position 


Tbol 


L boi 


Menv 


Class 


Bolocam ID 


230 GHz map? 


Other 




(SSTc2dJ...) 


(K) 


(L©) 


(M ) 








Names 


Ser-emb 1 


J182909.24+003132.3 


39 


4.1 


3.1 


Class 


Bolo 15 


Y 




Ser-emb 2 


J182952.44+003611.7 


42 


2.5 


1.3 


Class 


Bolo 24 


N 




Ser-emb 3 


J182854.84+002952.5 


51 


6.6 


3.8 


Class 


Bolo 8 


N 




Ser-emb 4 


J183000.72+011301.4 


54 


3.1 


6.6 


Class 


Bolo 28 


Y 




Ser-emb 5 


J182854.84+001832.6 


56 


0.4 


0.5 


Class 


Bolo 7 


Y 




Ser-emb 6 


J182949.56+011521.9 


56 


28.0 


20.0 


Class 


Bolo 23 


Y 


FIRS1 


Ser-emb 7 


J182854.12+002929.9 


58 


7.9 


4.3 


Class 


Bolo 8 


Y 




Ser-emb 8 


J182948.12+011644.9 


58 


5.4 a 


9.4 


Class 


Bolo 22 


Y 


S68N 


Ser-emb 9 


J182855.92+002944.7 


66 


2.5 


6.6 


Class 


Bolo 8 


N 




Ser-emb 11 


J182906.72+003034.3 


77 


4.8 


3.6 


Class I 


Bolo 14 


Y 




Ser-emb 15 


J182954.24+003601.3 


100 


0.4 


1.3 


Class I 


Bolo 24 


Y 




Ser-emb 17 


J182906.36+003043.2 


120 


3.8 


3.6 


Class I 


Bolo 14 


Y 





Note. — Position is the Spitzer source position. and L^i are calculated from the broadband SEP from 1.25 /im 

to 1100 fim (2MASS, IRAC, MIPS, Bolocam) and M env from the total 1.1 mm flu x, as described inlEnoch et all (2009a). 
Note that both L-^d and M env are larger by a factor of 2.5 than those listed in Enoch et al. (2009a) due to the larger 
cloud distance assumed here. "Bolocam ID" denotes the Bolocam 1.1 mm core with which the source is associated; note 
that more than one embedded protostar is often associated with a given Bolocam core, due to the 31" Bolocam resolution. 
"230 GHz map?" indicates if the source is included in the 230 GHz study. 
a Source is saturated at 70/im, so L^ \ is likely an underestimate. 



These studies have found that disks are well established 
in Class I, with typical masses of 0.01 — 0.5 M . Although 
circum-protostellar disks have been detected in a num- 
ber of Class sources, most detailed studies of Class 
sources have been limited to the brightest or most well- 
known objects due to instrumental limitations, a lack of 
complete target samples, and the difficulty of separating 
d isk and envelope emission in such young sources. 

|j0rgensen et al.l (j2007l l2QQ9h have recently completed 
a large Submillimeter Array (SMA) survey of ten Class 
and ten Class I sources, findin g similar disk masses 
in Class and Class I. Although U^rgensen et al.l ([2009) 
see no clear evolutionary signature in disk flux or mass, 
they do see an increase in the ratio of disk mass to enve- 
lope mass with evolutionary stage. The SMA study still 
relies on well known sources from different regions, how- 
ever, rather than a complete or well defined sample. A 
study characterizing "typical" Class disks, in a sample 
that is complete down to some well defined luminosity or 
envelope mass limit, is still needed. 

We have recently begun a program to characterize 
disks and envelopes in a complete sample of Class pro- 
tost ars in nearby low mass star- forming regions. Com- 
bining CARMA 230 GHz continuum imaging and Spitzer 
IRS mid-infrared (MIR) with radiative transfer model- 
ing of this sample will help to address several fundamen- 
tal questions about the structure and evolution of the 
youngest protostars: 1) How soon after the initial col- 
lapse of the parent core does a circum-protostellar disk 
form? 2) What fraction of the total circum-protostellar 
mass resides in the disk, and does this fraction vary with 
time? 3) Are large " holes" in the inner envelope (e.g. 
Ij0rgensen et a l. 2005) common at early times? 

Millimeter maps and MIR spectra provide complemen- 
tary approaches to these questions. The amount of flux 
escaping at A < 50 fim from embedded protostars is sen- 
sitive to the opacity clos e to the protostar, and thus the 
envelope structure (e.g. lJ0rgensen et al.l 12005). While 
the MIR flux is insensitive to disk properties, high res- 
olution millimeter continuum mapping directly detects 
emission from dust grains in the disk. Millimeter ob- 
servations with good ^-coverage and calibration, com- 



bined with radiative transfer models, can separate the 
disk from the envelope and constrain the disk mass and 
size. 

Our ultimate goal is to characterize the disk mass, size, 
and inner envelope structure of typical low-mass Class 
protostars, and to quantify any tren ds with evolutionary 
indicators. In lEnoch et al.l (j2009bl ). hereafter Paper I, 
we presented results for Serpens FIRS 1, a well known 
Class source, as a pilot for the full program. We found 
that disks can be quite massive at very early times; radia- 
tive t ransfer modeling (RADMC; Dull emond fc Dominikl 
2004) compared to CARMA visibilities and spectral in- 
formation indicated that FIRS 1 has a disk with mass 
Mdisk = 1.0 M and radius Rdisk — 300 AU, making up 
approximately 13% of the total circum-protostellar mass 
in a source that is about 10 5 yr old. 

Here we present high-resolution millimeter imaging of 
9 embedded protostars in the Serpens Molecular Cloud, 
including 6 of the 9 known Class sources in the cloud. 
The Combined Array for Research in Millimeter- Wave 
Astronomy (CARMA0) is an ideal instrument to study 
the disk properties of deeply embedded protostars. With 
15 antennas operating at A = 1 mm, CARMA provides 
the excellent ^-coverage required to separate disk emis- 
sion from envelope emission in young protostars. 

In Section [2] we describe the Serpens sample, and in 
Section [3] outline the CARMA observations and data re- 
duction. CARMA maps and visibilities are presented in 
Section HI followed by an initial estimate of disk masses 
from the visibility amplitude at large uv -distances in Sec- 
tion [5l Individual sources are discussed in Section [6l 
We comment on the observed multiplicity fraction in our 
sample in Section [3 

2. SERPENS SAMPLE 

Recent large surveys of nearby molecular clouds 
at mid-infrared and (sub) millimeter wavelengths have 
made it possible to define complete samples of Class 
protostar s based on lumino s ity or envelope mass 
limits (e.g. lHatchell et~aTI 120071 : Uorgensen et al.l 120081 : 

2 http://www.mmarray.org/ 
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Dunham et "all 120081 : Enoch et all l2QQ9al : lEvans et all 
2009D . Our study is based on the complete sample 
of 39 Class protostars in the Serpens, Perseus, and 
Ophiuc hus molecular clouds, identified by lEnoch et al.l 
(20093), identified by comparing large-scale Spitzer 
IRAC and MIPS with Bolocam (jGlenn et al.l I2003D 
1.1 mm continuum surveys of the three clouds. Here we 
focus on the Serpens cloud, which contains 34 embedded 
protostars, 9 Class sources and 25 Class I, and is com- 
plete to envelope masses M env > 0.25 M and internal 
luminosities L int > 0.05 L . 

The Serpens molecular cloud is an active star forming 
region with high extinction and a high density of young 
stellar objects (YSOs). It has been studied extensively at 
near-infrared, far- infrared, submillimeter, and millimeter 
wavelengths (e.g. Eiroa & Casali 1992; Hurt & Barsonv 
119961: lLarsson et a l. 2000; Davis et al. 1999; C asali et al.l 
119931 : iTesti fc Sargentl I1998D . We adopt a distance of 
d = 415 ± 25 pc based recent VLBA trigonometric par- 
allax measurements of a young stellar object in Ser- 
pens (|Dzib et al.ll2010f ). Note that this is significantly 
larger than the d = 260 pc that we have assumed in our 
previo us work (Paper I; lEnoch et al.l [20071 : lEnoch et al.l 
2009a). To compare to our previous disk and envelope 
masses in Serpens, masses given here should be decreased 
by a factor of 2.5 and physical sizes by 1.6. Mass ratios 
will remain unchanged. 

Targets a re listed in Tabled] in cluding the "Ser-emb#" 
name from lEnoch et al.l (|2009af ). Spitzer position, bolo- 
metric luminosity (Lb i), bolometric temperature (T& j), 
and total envelope mass (M env ). Our goal here is to im- 
age all known Class protostars in Serpe ns with high 
resolut ion and high fidelity with CARMA. lEnoch et al.l 
(20093) found 9 Class protostars in Serpens. Three of 
the 9 known Class source were not included in the 
230 GHz observations, however. These three sources 
were either not detected, or only marginally detected, 
in preliminary 110 GHz lower-resolution CARMA maps 
(see Section [A]). 

All undetected sources are in confused regions, with 
multiple Spitzer sources associated with a single Bolo- 
cam 1.1 mm core; at least one source in each core is 
detected and included here. The undetected sources do 
not have particularly low luminosity or low envelope mass 
(see Table [T]), so excluding them does not bias the survey 
by luminosity or circum-protostellar mass. In addition to 
six Class targets, we include three Class I sources that 
are near the Class 0/1 cutoff. 

Two sources in our s ample (Ser-emb 4 and em b 6) were 
observed as part of the Hogerheijd e et all (jl999[ ) interfer- 
ometric study of four deeply embedded protostars in Ser- 
pens, but of the two only Se r -emb 6 was detected. The 
other two iHogerheijde et al.l (|1999f ) sources were classi- 
fied as Class I by Spitzer. 

As described in lEnoch et al.l (j2009af ). L^oi and T^i are 
calculated from the spectral energy distribution (SED) 
from 1.25 jam to 1100 jam. Broadband infrared data 
for our targets are taken from the "From Molecular 
Cloud Cores to Planet Forming Dis ks" Spitzer Legac y 
program ("Cores to Disks" or c2d; lEvans et al.l 12003). 
which image d one square degreee in the cloud with IRAC 
and MIPS (jHarvev et al.l I2006L 1200/1 ). The same re- 
gion was mapped at A = 1.1 mm with the Bolocam 
bolometer array at the Caltech Submillimeter Observa- 
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Fig. 1. — Overview of the Serpens Main cluster and Cluster B 
region, with a finding chart for embedded protostars observed with 
CARMA. The Spitzer 24 /xm map from c2d is shown in grayscale, 
with Bolocam 1.1 mm continuum contours overlaid (contours are 
50, 100, 200, 400, 600. ..1400 mJy beam" 1 ). Although the Serpens 
cloud, as defined by the Ay = 6 contour, extends beyond the 
region shown, all embedded protost ars are fo und within this area. 
"emb#" designations are from Enoch et al. (2009a|), and sources 
marked with a "(*)" are Class protostars not included in the 
230 GHz study (see Sections [2] and [Aj. 
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Fig. 2. — Spitzer three color (8 jum, 24 /im, 70 /im) images of the 9 embedded protostars imaged here with CARMA. Bolocam 1.1 mm 
continuum contours are shown at 50,100,150,200,400,600,800,1000 mJy beam -1 . Each target is centered; Ser-emb 11 and emb 17 are 
< 15 /; apart and are shown in the same panel. 



tory (|Enoch et al.ll2007[ ). Together with public 2MASS 
catalogs, these data provide wavelength coverage at 1.25, 
1.65, 2.17 /im (2MASS), 3.6, 4.5, 5.8, 8.0 /im (IRAC), 24, 
70, 160 /im (MIPS), and 1100 /im (Bolocam). 

The envelope mass M env is calculated from the to- 
tal A = 1.1 mm single dish flux, assuming the en- 
velope is optically thin at 1.1 mm, a dust opacity of 
Kxmm = 0.0114 cm 2 g _1 (|Ossenkopf fc Hennnigll 19941 ). a 
dust temperature of Tjj = 15 K, and d = 415 pc. Here we 
treat the envelope mass as the total circum-protostellar 



mass, as it may include some flux contribution from a 
disk component. 

Figure [1] gives an overview of the Serpens environment 
with Spitzer 24 /im (grayscale) and Bolocam 1.1 mm 
continuum (contours) maps. The locations of our embed- 
ded protostar targets are indicated. Most are associated 
wi th the Serpens Main cluster to the nor th or Cluster 
B (jDiupvik et al,ll2QQ6t lHarvev et aLll2006D to the south, 
but a few are more isolated. Three color Spitzer images 
(8, 24, 70 /im) of the 9 targets imaged with CARMA are 
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shown in Figure [2j with Bolocam 1.1 mm contours as in 
Figure [U SEDs of each source from 1.25 /im to 1100 jam 
can be found in lEnoch et al.l (|2009a[ ). 



3. OBSERVATIONS AND DATA REDUCTION 

Continuum observations at v = 230 GHz (A = 1.3 mm) 
were completed with CARMA, a 23 element interfer- 
ometer consisting of nine 6.1-m, six 10.4-m, and eight 
3.5-m antennas. The 6.1-m and 10.4-m antennas were 
used to obtain 230 GHz continuum observations in the B 
(100-1000 m baselines), C (30-350 m), D (11-150 m), 
and E (8 — 66 m) antenna configurations between Octo- 
ber 2007 and January 2010. Three to four sources were 
observed in each track; although the ^-coverage is not 
identical for each target, every source was observed in 
each configuration, and at least ±3 hours around tran- 
sit. Integration times varied based on the expected flux 
from single dish observations. 

Three correlator bands were configured for continuum 
observations with 468.75 MHz bandwidth, for a total 
bandwidth of 2.82 GHz (two sidebands per band). A 
bright quasar (1751+096) was observed approximately 
every 15 — 20 minutes, to be used for complex gain cali- 
bration. For B configuration tracks, a secondary calibra- 
tor (1830+063) was also observed every 15 minutes. Ab- 
solute flux calibration was accomplished using 5 minute 
observations of Uranus, Neptune, or MWC 349. The 
overall calibration uncertainty is approximately ±20%, 
from the reproducibility of the phase calibrator flux on 
nearby days. A passband calibrator, typically 3C454.3, 
was observed for 15 minut es during each set of obser- 
vations, and either optical (jCorder et al.l [2010) or radio 
pointing was performed every one to three hours. 

Data from different array configurations were com- 
bined to provide excellent iw-coverage from 4 to 500 /cA, 
providing sensitivity to physical scales from 170 AU to 
2 x 10 4 AU. Small 7-point mosaics were made in the com- 
pact configurations (D and E) in order to mitigate spa- 
tial filtering by the interferometer. Mosaicing provides 
better sensitivity to the spatially extended protostellar 
envelopes. 

Observations in the most extended B configuration uti- 
lized the CARMA Paired Antenna Calibration System 
(C-PACS) to correct for p hase variations on sub-minute 
timescales ([Perez et al.ll2010f ) . For the longest B configu- 
ration baselines, C-PACS pairs 6.1-m and 10.4-m anten- 
nas (science array) with 3.5-m antennas (reference array) 
that operate at A = 1 cm. The reference array continu- 
ally monitors the atmospheric phase using a strong point 
source close to the science target (here 1751+096, also 
the gain calibrator, which is ~ 13deg away), while the 
science array slews between the science targets and gain 
calibrator. The delay measured by the reference array 
is later used t o correct the scien ce data for atmospheric 
variations (see lPerez et al.ll2010l for details). The 3.5-m 
antennas were tuned to 30 GHz, with a 14-band wide- 
band correlator providing 7 GHz of bandwidth. 

Calibration and imaging are accomplished 
with the MIRIAD data reduction package 
(jSault. Teuben. fc Wrightl 119951 ). Calibration of 
each track consists of applying line-length calibration, 
passband correction, phase correction on 15 — 30 min 



intervals, amplitude calibration on 60 min intervals, and 
absolute flux calibration. 

For B configuration data, delays are calculated from 
the 30 GHz observations of 1751+096 and applied to each 
paired baseline on 4 s intervals after passband calibra- 
tion. The 230 GHz gain calibrator observations are then 
used to apply longer- interval (15 min) phase corrections, 
to remove slowly varying instrumental delay differences 
between the reference and science arrays. C-PACS cor- 
rections usually work best for close (< 6 deg) atmospheric 
calibrators (Zauderer et al., in preparation); given that 
1751+096 is 13 deg from our science targets, we only ap- 
ply C-PACS calibration for tracks where it significantly 
improves phase coherence on the secondary calibrator 
(1830+063, 5 deg from the science targets). When ap- 
plied, C-PACS calibration improves the signal to noise 
of the peak flux of 1830+063 by a factor of 1.5 — 6. 

Images are formed by inverting the calibrated visibili- 
ties using multi-frequency synthesis with natural weight- 
ing to optimize sign al to noise, cle aned with a Steer 
CLEAN algorithm (jSteer et al.lll984fl . and restored with 
a Gaussian fit to the synthesized beam. Beam sizes for 
the combined maps (including all configurations) are ap- 
proximately 1", corresponding to a physical resolution of 
415 AU. 

We also calculate the azimuthally averaged visibility 
amplitude as a function of ^-distance. Due to mosaic- 
ing in the compact arrays, visibilities must be combined 
over multiple pointing centers, and are calculated in three 
steps. The visibilities are first modified in phase to ac- 
count for any offset between the source position and the 
pointing position. The visibilities are then divided by 
the primary beam response at the location of the peak 
of the flux in the various pointings. Finally, we plot the 
weighted mean of the visibilities over the different point- 
ing positions. With this method, very extended compo- 
nents or bright multiple sources can introduce errors in 
the visibility amplitudes. However, if the emission is very 
compact relative to the 30" primary beam, the error in- 
troduced by this scaling of the visibilities from the offset 
pointings is minimal. 

4. MAPS AND VISIBILITIES 

The resulting 230 GHz maps are shown in Figure [3] 
In addition to maps including all data, two maps includ- 
ing only m;-distances < 50 kA and > 50 kA are shown, 
to better illustrate any extended and compact emission, 
respectively. 

The 230 GHz peak flux, total flux, and deconvolved 
size of each detected source, determined from a Gaussian 
fit (MIRIAD's imfit) to the combined map, are given 
in Table El The rms noise level in the central region, 
excluding sources, is also listed (ranging between 0.9 — 
8 mJy beam -1 ), along with the synthesized beam size, 
or resolution (typically ~ 1" or 415 AU). Multiple entries 
are given when more than one source is detected in the 
CARMA map. 

Plots of the visibility amplitude versus ^-distance for 
each embedded protostar are shown in Figure [6l Ex- 
tended envelope structures, as expected for embedded 
Class protostars, are visible as an amplitude peak at 
m;-distances < 20 kX in the visibility plots (Figure [6]). 
Only Ser-emb 15 does not show evidence for a strong 
peak at ^-distance < 20 kA. The peak visibility am- 
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TABLE 2 

CARMA 230 GHz source properties 



ID 


230 GHz position 


Peak flux 


Total flux 


rms 


beam 


deconvolved size 






(mjy bin -1 


) (mJy) 


(mJy bin - 1 ) 


(arcsec) 


(arcsec) 


Ser-emb 1 


18:29:09.09 +00:31:30.9 


94.7 


215 


3.0 


1.2x1.1 


1.5x1.1 


Ser-emb 4 (N) 


18:30:00.30 +01:12:59.4 


6.3 


245 


1.5 


1.2x1.1 


8.2x6.3 


Ser-emb 4 (E) 


18:30:00.70 +01:12:56.2 


6.1 


30.3 


1.5 


1.2x1.1 


2.6x2.1 


Ser-emb 4 (S) 


18:30:00.62 +01:12:49.5 


5.6 


130 


1.5 


1.2x1.1 


7.2x4.3 


Ser-emb 5 


18:28:54.91 +00:18:32.4 


6.1 


14.3 


0.8 


1.3x1.2 


1.6x1.2 


Ser-emb 6 


18:29:49.79 +01:15:20.4 


423 


1390 


6.7 


0.9x0.9 


1.6x1.2 


Ser-emb 7 


18:28:54.04 +00:29:29.7 


9.7 


281 


1.9 


1.0x1.0 


5.9x4.6 


Ser-emb 8 


18:29:48.07 +01:16:43.7 


40.1 


3970 


4.3 


0.7x0.6 


7.0x5.9 


Ser-emb 11 (W) 


18:29:06.61 +00:30:34.0 


25.6 


107 


2.8 


1.0x1.0 


1.9x1.6 


Ser-emb 11 (E) 


18:29:06.75 +00:30:34.4 


13.3 


55.2 


2.8 


1.0x1.0 


1.8x1.6 


Ser-emb 15 


18:29:54.30 +00:36:00.8 


41.5 


77.5 


1.7 


1.0x1.0 


1.1x0.8 


Ser-emb 17 


18:29:06.20 +00:30:43.1 


42.2 


109 


2.8 


1.0x1.0 


1.3x1.2 



Note. — Peak fluxes, total fluxes, and source sizes are determined by a Gaussian fit in the combined 230 GHz maps. 
The 1 a rms is cacluated in the central 15" of each map, excluding bright sources. An angular size of 1" corresponds to a 
physical size of 415 AU. 



plitude is equivalent to the total 230 GHz flux of the 
sourc^E The interferometer does filter out flux at uv- 
distances less than 4 kX (physical scales > 2 x 10 4 AU), 
corresponding to the separation of the closest antenna 
pairs, so the total flux observed by CARMA will often 
be less than the total single dish flux. 

A point source, such as an unresolved disk, will ap- 
pear as a constant amplitude in the visibility plot, while 
a resolved disk will contribute an additional wide Gaus- 
sian component to the visibilities. Although the most 
common interpretation of a compact visibility compo- 
nent is a disk, without dynamical information or inde- 
pendent constraints on the envelope structure, compact 
emiss ion could also be a ttributed to a power law enve- 
lope (Loonev et al. 2003) or non-axisymmetric envelope 
structure (jLooney et al.l 12007). With these caveats in 
mind, here we interpret compact emission as evidence of 
a circum-protostellar disk. 

We can see immediate evidence for a compact disk 
component, as non-zero amplitudes at ^-distances > 
50 &A, in most sources (e.g. Ser-emb 4, emb 7,emb 11, 
emb 17, emb 15). Other sources have significant, but 
non-constant, flux for uv- distance 50 — 100 /cA, which 
may be due to a partially resolved disk (e.g. Ser-emb 1, 
emb 6, emb 8). 

5. ESTIMATED DISK MASSES 

Accurately measuring disk masses requires compar- 
ing the 230 GHz visibilities and spectral informa- 
tion with detailed radiative transfer models, as in 
Paper I. Radiative transfer models such as RADMC 
(Dul lemond fc Dominikll2004f ) simultaneously model the 
envelope and the disk, taking into account any enve- 
lope component at small scales, as well as a possibly 
resolved disk. In addition, radiative transfer models self- 
consistently solve for the temperature as a function of po- 
sition. Work to constrain the disk and envelope structure 
of these sources with grids of radiative transfer models is 
currently ongoing (Enoch et al., in prep.). 

For now, we make a first estimate of the disk mass 
using the 230 GHz flux at a ^-distance of 50 kA. An 

3 Note that this can be larger than the total flux listed in Table [2] 
because a Gaussian fit does not always capture the full extent of 
an extended source. 



unresolved disk will have a constant flux independent 
of ^-distance, while the resolved envelope flux falls 
off quickly with distance. Radiative transfer mod- 
els of interferometric observations of a number of em- 
bedded protostars have indicated that protostellar en- 
velopes con tribute little emission at baselines longer than 
50 kA (e.g. iKeene fc Masson I Il990l: lHarvev et all 120031 : 
|j0rgensen et al.ll2005l ). Using similar models to separate 
the disk and envelope in Ser-emb 6, we find that the en- 
velope flux for ^-distances > 30 kA is nearly zero (see 
Figure 6 of Paper I), contributing < 2% of the total flux 
at 50 kA. 

The choice of 50 kA, as was adopted bv lJ0rgensen et aTl 
(2009), is a good compromise between minimizing the 
envelope contribution and maximizing the signal-to-noise 
of the disk, which falls off with increasing baseline length 
(due to decreased uv coverage). 

To calculate the mass, we assume that the disk is op- 
tically thin at v = 230 GHz: 



Mdisk 



d 2 S> 



50k\ 



^1.3ram(^D)ftl.3ri 



(i) 



where S$o k\ is the flux at 50 kA, d = 415 pc is the 
cloud distance, and £?i.3 mm is the Planck function for 
a du s t temperature of Tp = 30 K (e.g. IVisser et al.l 
2009; lJ0rgensen et al.ll2009l ). We assume the dust com- 
position is similar in the disk and envelope, taking 
the dust opacity per gram of gas at A = 1.3 mm, 
fti.smm = 0.009 cm 2 g -1 , from Table 1 column 6 of 
lOssenkopf fc Henningj (| 1994ft , for dust grains with thin 
ice mantles. 

Here we assume that there is a negligible contribution 
by the envelope to the flux at 50 kA. If approximately 
30% of the flux a t 50 kA can be co ntributed to the enve- 
lope, as found by J0r gensen et al.l (|2QQ9h based on radia- 
tive transfer models of 8000 AU envelopes, disk masses 
could decrease by 30%. 

On the other hand, if any of the disks are partially 
resolved (i.e. have radii larger than 170 AU), the flux at 
50 kA will underestimate the disk mass. For Ser-emb 6, 
the flux at 50 kA suggests a mass that is 30% smaller 
than the mass from detailed radiative transfer modeling 
([Enoch et al.l l2QQ9bh . Furthermore, if any part of the 
disk is optically thick, or the outer regions of the disk 



Disks in Class Protostars 



7 



i/IO 




Ser-emb 1 
Center: 18 29 09.0 +00 31 31.£ 




PI// . 



' 



10 -10 -20 
Offset (") 



750 AU 7" ( b ) 

_l I I I I I £_l 2_l I I I I I I L 



5 0-5 
Offset (") 




-10 



4 3 2 1 -1-2-3-4-5 
Offset (") 



CD 
if) 

O 








300 AU 



(c) 

■ /T^s 



20 



10 -10 -20 

Offset (") 



1 o 












: % 




a i— i o 

1500 ALL 






10 -10 -20 

Offset (") 



10 


5 0-5 
Offset (") 


-10 -1 


, , , ''i i t i .;■ f 

Ser-emb 5 
'. Center: 18 28 54.8 +00 18^5Z^ 




















m 1 — 1 

. • 750 AU 

, , , , i , , , , i , , , , 






5 
Offset (") 


-5 -1 



-2-3-4-5-6-7-8-9 
Offset (") 




10 



4 3 2 1 -1-2-3-4-5 
Offset (") 



Fig. 3. — CARMA 230 GHz images of Serpens deeply embedded protostars. The combined maps, shown in panel (b), include data from 
B, C, D, and E configurations, with 7-point mosaics in D and E. Large-scale and small-scale emission maps are shown in panels (a) and 
(c), respectively; maps in panel (a) are made including only i//u-distances < 50 kA, and maps in panel (c) with u^-distances > 50 kA (with 
the exception of Ser-emb 4, which includes ^-distances > 30 kA due to low signal-to-noise on long baselines). Contours are shown at 
3,5, ..15, 20. ..60, 80. ..140 cr, where a is given in Table [2] with lighter contours starting at 15 a and negative contours indicated by dotted 
lines. Synthesized beam sizes (approximately \" or 400 AU for combined maps) are indicated by solid ellipses. Offsets from the pointing 
centers are likely due to errors in the Spitzer positions, which are determined by the shorter wavelength data where outflows can cause 
extended or offset emission. 
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Fig. 4. — Figure 3 cont. 
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Fig. 6. — Azimuthally averaged visibility amplitudes as a function of m;-distance for the Serpens embedded protostars. The value 
expected for zero signal, or bias, is shown for reference (dotted line). Averaged visibilities are vector averaged at the measured source 
position (Table [2]), and not the pointing center. Error bars reflect the variable m;- cover age, and atmospheric decorrelation on the longest 
baselines (> 100 or > 200 kA depending on source flux); they do not include the ~ 20% absolute calibration uncertainty. Visibilities are 
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TABLE 3 
Estimated disk masses 



ID 


Flux(50 kA) 


Estimated M^isk 




(mJy) 


(M ) 


Ser-emb 1 


120 


0.28 


Ser-emb 2 




< 0.02 


Ser-emb 3 




< 0.02 


Ser-emb 4 (N) 


26 


0.05 


Ser-emb 5 


17 


0.04 


Ser-emb 6 


750 


1.7 


Ser-emb 7 


66 


0.15 


Ser-emb 8 


110 


0.25 


Ser-emb 9 




< 0.02 


Ser-emb 11 (W) 


52 


0.13 


Ser-emb 15 


64 


0.15 


Ser-emb 17 


69 


0.15 



Note. — The flux at 50 kA is the average flux for 
^-distances 40 — 60 kA. Disk masses are estimated from 
Equation [TJ and upper limits for Ser-emb 2, emb 3, and 
emb 9 are based on preliminary 110 GHz maps (Sec- 
tion [Aj. Uncertainties in disk masses are at least 50%, 
as discussed in the text. 

have a dust temperature well below 30 K, the true mass of 
the disk could be considerably higher than our estimates. 
We therefore assign uncertainties to our estimated disk 
masses of at least ±50%. 

Estimated disk masses are given in Table [3j and are 
plotted as a function of bolometric temperature, a mea- 
sure of the evolutionary state, in Figure [9j We do not see 
a trend with evolutionary state in the narrow T^i range 
observed. The mean disk mass in our 9 source sample is 
0.33 M . 

We must also account for the three Class sources 
not included here due to non-detections in preliminary 
110 GHz maps. If Ser-emb 2, emb 3, or emb 9 had a 
compact disk with Mdisk ^ 0-02 M they would have 
been detected at 3mm (see Section [A]), so we assume that 
these sources do not have significant disks. Thus, we find 
evidence for circum-protostellar disks in 6 of the 9 Class 
sources in Serpens. If we include these non-detections, 
we find a mean disk mass < Mdisk >= 0.2 M . 

The disk to envelope mass ratio (Mdisk/ M env ) is also 
shown versus T^i in Figure [9j where M env is from Ta- 
ble [2j Disk masses are between 1% and 10% as large as 
envelope masses in our sample. The mean Mdisk / M env is 
6%. Again we see no clear trend in the disk to envelope 
mass ratio with T^i within this young sample. 

Our high detection rate in this complete sample sug- 
gests that disks must form very early in the star for- 
mation process. The Class phase l ikely lasts approx- 
imately 0.2 Myr (En och et al.ll200 9a). and we have de- 
tected disks in 6 of the 9 Class targets, which repre- 
sent all Class sources with L^ nt > O.O5L and M env > 
0.25M env in the approximately 3.5 deg cloud area. Thus 
our results suggest that 2/3 of embedded protostars in 
Serpens have formed an Mdisk > 0.02 M disk within 
the first 0.2 Myr. 

Although there is no overlap between the two sam- 
ples, the disk masses for our co mplete Serpens s ample 
are roughly consistent with the J0rg ensen et al.l {2009) 
Class disk masses. We also see similar disk to envelope 
mass ratios in Class 0, ranging from 1 — 10% in both 
samples. 



|j0rgensen et al.l ([2009) found average disk masses 
in bo th Class and Clas s I of Mdisk = 0.05 M , 
while lAndrews fc Williamsl ([2007) found a lower aver- 
age Class I disk mass in Ophiuchus (0.015 M ). Our 
mean Class disk mass is somewhat higher, 0.2 M , but 
this may be a result of temperatur e variations with lu- 
minosity (e.g. J0rgen sen et al.ll2009[ ), or the fact that we 
have not corrected for the envelope contribution at 50 kA. 
Given the uncertainties, ou r results appear t o be c onsis- 
tent with the conclusion of J0rge nsen et al.l ([2009), that 
disk masses are fairly constant through the embedded 
phases of protostellar evolution. 

Such a scenario i s supported by r ecent hydrodynam- 
ics simulations by IVorobvovl ([20091 ), which find time- 
averaged disk masses in Class and Class I of ~ 0.1 M 
for both viscous and self-gravitational disks. More ro- 
bust measurements of disk masses and similar studies of 
the Class populations in other clouds are needed to 
confirm these results. 
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Fig. 9. — Estimated disk masses (top) and ratio of the estimated 
disk mass to envelope mass (bottom), as a function of bolometric 
temperature. Upper limits for three Class sources in our sample 
are based on preliminary 110 GHz maps (Section EJ- No strong 
trend with is seen for either the disk mass or disk to envelope 
mass ratio. 

6. INDIVIDUAL SOURCES 

6.1. Ser-emb 1 

Ser-emb 1 is located in Serpens Cluster B (Figured]), 
and is associated with the sing l e dish Bolocam 1.1 mm 
core Ser-Bolo 15 (jEnoch et al.l f2QQ7h . Ser-emb 1 wa s 
identified as a Class source by (jDjupvik d~aTl 2006). 
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Ser-emb 1 has the lowest bolometric temperature (39 K), 
and thus presumably is the least evolved, of all the proto- 
stars in Serpens. It is relatively isolated; a single Spitzer 
mid-infrared source is associated with the approximately 
10 4 AU diameter Bolocam core, and the nearest proto- 
star is - 70" (3 x 10 4 AU) away. 

The 230 GHz map of Ser-emb 1 (Figure [3]) shows a sin- 
gle compact source, which is centered on the mid-infrared 
source. The total flux from a Gaussian fit is 215 mJy (Ta- 
ble EJ), but the peak in the visibility plot (Figure [6j right) 
indicates a total flux in the CARMA data of 560 mJy. 
After scaling the 270 GHz (A = 1.1 mm) Bolocam flux 
to 230 GHz, assuming a dust spectral index of j3 = 1.5, 
we find that CARMA recovers approximately 95% of the 
single dish flux. 

The visibility plot indicates substantial flux at inter- 
mediate ^-distances (30 — 100 kA), which could be ex- 
plained by either a resolved disk or a very compact en- 
velope. The disk mass estimated from the flux at 50 kA 
is 0.28 Mq, or 9% of the total circum-protostellar mass, 
although this may be severely underestimated if the disk 
is partially resolved. 

6.2. Ser-emb 4 

Ser-emb 4 is located in the Serpens Main cluster, in a 
crowded region (several other protostars within 1'). It is 
associated with the Bolocam core Ser-Bolo 28, as are Ser- 
emb 19 and Ser-emb 29. Most of the single dish flux, and 
thus most of the envelope mass, is likely associated with 
Ser-emb 4, which has a much lower T^i (54 K compared 
to 129 K and 374 K, respectively). 

The 230 GHz map reveals three distinct sources; the 
northernmost is associated with the Ser-emb 4 mid- 
infrared source, while the other two appear to have no 
counterpart at A < 70 fim. The eastern and southern 
sources are located 10" (4200 AU) and 15" (6200 AU) 
from the primary northern source, respectively. 

Ser-emb 4 (N) has a diffuse extended envelope that 
is partially resolved out by the small synthesized beam. 
The extended nature of the envelope is also apparent in 
the narrow peak in the visibilities at small distances. 
The peak visibility amplitude is 250 mJy, only 19% of 
the total single dish flux. This difference is accounted 
for by the fact that the Bolocam core breaks into several 
sources with the higher resolution of CARMA. A con- 
stant amplitude of ~ 25 mJy is seen out to ^-distance 
~ 100 kA, suggesting an unresolved disk component with 
Mdisk ~ 0.05 M , approximately 1% of the total circum- 
protostellar mass. The two additional sources in the 
CARMA map also likely contribute to the single dish 
flux, making this 1% fraction a lower limit. 

6.3. Ser-emb 5 

Ser-emb 5 is located south of Cluster B, and is one 
of the few embedded protostars not in one of the two 
main clusters. It is associated with the extended Bolo- 
cam core Ser-Bolo 7, as is Ser-emb 25, but the two sources 
are separated by > V '. Ser-emb 5 has a low bolomet- 
ric luminosity (L^i = 0.4 L ) and low envelope mass 
(M env = 0.6 M ). It is the only Class source in Ser- 
pens that is considered a Very L ow Luminosity Object 
(VeLLO; see iDunham et aLl l2QQ8). VeLLOs are believed 
to be either very low mass protostars (i.e. proto-brown 



dwarfs), or stella r mass protostars w ith very low mass 
accretion rates ([Dunham et al.ll2008f ). 

The 230 GHz map contains a single, relatively com- 
pact source. The peak visibility amplitude is 68 mJy, 
accounting for 60% of the single dish flux. The flux at 
50 kA suggests a disk mass of 0.04 M , or 6% of the total 
circum-protostellar mass. Ser-emb 5 has the lowest disk 
mass of our nine 230 GHz targets. 

6.4. Ser-emb 6 

Ser-emb 6 is a well k nown Class protos tar also known 
as Serp ens FIRS 1 (lHaryev et al.l Il984f ) and Serpens 
SMM 1 (Cas ali et al.lll993[ ). It is located in the Serpens 
Main cluster and is associated with the Bolocam core 
Ser-Bolo 23. Ser-emb 6 is the brightest Class protostar 
in Serpens (L^i = 28 L ), and has the largest envelope 
mass (M env = 20 M ) of any embedded protostar in the 
region. 

The 230 GHz map of Ser-emb 6 shows a complex 
source with both compact and extended emission. Ap- 
proximately 58% of the single dish flux is recovered by 
CARMA, with a peak visibility amplitude of 2.3 Jy. In 
Paper I we compared the CARMA visibilities, Spitzer 
IRS spectra, and broadband SED to a grid of de- 
tailed radiative transfer mode ls calculated with RADMC 
(jDullemond fc Dominikll2004f ). The visibilities were best 
fit by a massive disk (Mdisk = 2.5 M when scaled for 
d = 415 pc), making up 13% of the total circumstellar 
mass and accounting for most of the flux at ^-distances 
30 — 150 kA. Using only the flux at 50 kA results in a 
disk mass of 1.7 M , or 9% of the circumstellar mass, an 
u nderestimat e of 30% compared to detailed modeling. 

iChoil ([2009) suggests that Ser-emb 6 is a proto-binary 
source with a separation of 800 AU. This second source 
corresponds to the small peak seen to the north-west 
of Ser-emb 6 (Figure [3]). Given its location along the 
outflow axis, and the clumpy nature of the envelope, we 
concluded in Pap er I t h at thi s emission is associated with 
the outflow, but Choi (j2009[ ) argues that the source is a 
second disk component based on the spectral index from 
3.5 cm to 7 mm. 

6.5. Ser-emb 7 

Ser-emb 7 is located in Cluster B, and is associated 
with the Ser-Bolo 8 core. Two other nearby protostars, 
Ser-emb 3 (30" away) and Ser-emb 9 (26" away) are asso- 
ciated with the same single dish core, but both were only 
marginally detected in initial 110 GHz CARMA map 
(Section [A]), suggesting that most of the core mass is 
associated with Ser-emb 7. Another possibility is that 
Ser-emb 3 and emb 9 have envelopes with a flat central 
density profile, which is resolved out by the interferom- 
eter. Neither Ser-emb 3 or emb 9 are included in our 
230 GHz study due to the lack of observed millimeter 
flux. 

The 230 GHz map of Ser-emb 7 shows a single ex- 
tended source that appears elongated in the north-south 
direction. Only 29% of the single dish flux is recovered 
by CARMA. The flux at 50 kA suggests a disk mass of 
0.15 M , or 4% of the total circum-protostellar mass, but 
the two nearby protostars may contribute to the single 
dish flux if they have envelopes that are resolved out by 
the interferometer. 
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6.6. Ser-emb 8 

Ser-emb 8, also known as S68 N (jMcMullin et al.lll994f ) 
or SMM 9 (jCasali et al.lll993l ). is located in the Ser- 
pens Main cluster. Although two Class I protostars 
are located within the FWHM of the Bolocam core Ser- 
Bolo 22, the core is centered on Ser-emb 8 and the major- 
ity of envelope mass appears to be associated with this 
object. 

The 230 GHz map indicates a compact source at the 
position of Ser-emb 8, with two additional sources ~ 12" 
and ~ 20" away. Approximately 20% of the single dish 
flux is recovered with CARMA; as with Ser-emb 4 and 
emb 7, nearby sources were likely blended in the Bolocam 
beam. The flux at 50 kA suggests a disk mass of 0.25 M , 
or 3% of the total circumstellar mass, likely a lower limit 
as nearby sources are contributing to the single dish flux. 

6.7. Ser-emb 11 and Ser-emb 17 

Ser-emb 11 and emb 17 are both Class I sources asso- 
ciated with the Ser-Bolo 14 core, located in Cluster B. 
The two protostars are separated by 12", or 5000 AU, 
suggesting that they may be physically associated. They 
have similar SEDs, but Ser-emb 17 is brighter at IRAC 
wavelengths, resulting in a highe r T^i (117 K compared 
to 77 K). In lEnoch et all (j2009afl . the total 1.1 mm flux 
of Ser-Bolo 14 was divided evenly between them to de- 
termine the envelope masses. 

The 230 GHz map reveals two compact sources at the 
positions of Ser-emb 11 and emb 17. Ser-emb 11 is fur- 
ther resolved into two components separated by 2.1" or 
870 AU, with a flux ratio of approximately 2. This is the 
only source in our sample with clear evidence for a binary 
or multiple component for separations < 2000 AU. 

The visibility amplitudes of both Ser-emb 11 and 
emb 17 have narrow peaks, indicating very extended en- 
velopes, and nearly constant amplitudes to ^-distance 
> 100 kA, suggesting compact disk components. To- 
gether, the peak visibility amplitudes of Ser-emb 11 and 
emb 17 account for approximately 60% of the total single 
dish flux. The fluxes at 50 kA are similar, resulting in disk 
masses of 0.13 M (emb 11) and 0.15 M (emb 17), to- 
gether approximately 4% of the total circumstellar mass. 
Note that visibilities and disk mass estimates are calcu- 
lated only for the brighter emb 11 component (W). 

6.8. Ser-emb 15 

Ser-emb 15 is a Class I source located to the east of 
Cluster B, in the Ser-Bolo 24 core. The Bolocam core is 
also associated with the Class source Ser-emb 2, which 
is 28" or 12000 AU away from emb 15. Ser-emb 2 was 
not detected in preliminary 110 GHz CARMA imaging 
(Figure [T0|) . As there is significant 1.1 mm emission asso- 
ciated with emb 2 (see also Figure [2]), the non-detection 
may be due to the lack of a compact disk combined with 
a flattened envelope density profile. 

The 230 GHz map of Ser-emb 15 shows a single com- 
pact source. The visibility amplitudes are nearly con- 
stant for ^-distances 5 — 200 kA, without the peak at 
small uv distances expected for an extended envelope. 
Rather, the visibilities are consistent with a bright unre- 
solved disk. The CARMA observations recover 54% of 
the total single dish flux associated with Ser-emb 15. The 
flux at 50 kA suggests a disk mass of 0.15 M , only 13% 



of the total circum-protostellar mass. If a large fraction 
of the 1 mm flux is coming from the disk, however, the 
envelope mass calculated assuming Tjj = 15 K is likely 
overestimated. Thus, Ser-emb 15 appears to have a com- 
pact disk surrounded by a relatively low mass envelope, 
as expected if Class I source s have exhau sted at least half 
of their envelope mass (e.g. Andre 1994). 

7. MULTIPLICITY 

The very high resolution millimeter maps presented 
here are also sensitive to binary and multiple sources, as- 
suming that each source has a separate disk or envelope 
component. The multiplicity fraction at very early times 
is essential for understanding the formation of multiple 
systems, but very little is currently known about multi- 
plicity in the Class stage. 

Although our observations do not have uniform sensi- 
tivity, the sensitivity to small scales (down to 250 AU) 
makes this sample a valuable addition to multiplicity 
studies. In only one source, Ser-emb 11, do we find strong 
evidence for multiplicity on scales < 2000 AU (5"). Ser- 
emb 11 (E) and (W) are separated by 870 AU. Although 
Ser-emb 11 is a Class I source, it is very near the Class 
cutoff (Tjjoi = 77 K), and is probably young enough to 
be included with the Class sample. For our sample 
of 9 sources, this yields a multiplicity ratio of approx- 
imately 10%. Ser-emb 6 may also be a binary source 
wi th 800 AU separation based on the VLA observations 
of iChoil (j2QQ9h . so it is possible that the fraction is as 
high as 20%. 

Approximately 21% of Class I protostars 
have companions in the 250 — 200 AU range 
(jConnellev. Reipurth. fc Tokunagal 1200 8). with sim- 
ilar binary frequencies in more evolved Class II (~ 16%; 
Kraus et al., in pr eparation) and main s e quence solar- 
type s tars (~ 10%: iRaghavan et al.ll20l(l . IMaurv et all 
(j201Q[ ) recently observed 5 Class protostars with 
sub-arcsecond resolution with IRAM-PdBI to study 
the fraction of binary and multiple systems down to 
50 AU separations. No Class binaries were found 
with separations < 1500 AU, in contrast to the ~ 20% 
multiplicity fractions observed in m ore evolved objects. 

Combining our sample with the IMaurv et al.l (j201Qh 
study, the resulting low Class multiplicity fraction 
(~ 10%) is gene rally consistent the suggestion of 
IMaurv et al.l (j2010f ). that the binary fraction in the 75- 
2000 AU range increases from the Class to Class I stage, 
possibly due to dynamical evolution of initially very tight 
(< 75 AU) or very wide (> 2000 AU) systems. With still 
relatively few sources, however, the statistical uncertain- 
ties are too large to conclusively determine that Class 
sources have fewer intermediate scale binaries than more 
evolved protostars. 

8. SUMMARY 

We have completed high-resolution, high-fidelity mil- 
limeter imaging of nine embedded protostars in the Ser- 
pens molecular cloud with the CARMA interferometer. 
Our sample is based on the complete (M env > 0.25 M , 
Lbol ^ 0.0 5 Lfp) census of em bedded protostars in Ser- 
pens from lEnoch et all (j2009af h and includes 6 of the 9 
Class sources in Serpens. 

230 GHz continuum observations were completed in B, 
C, D, and E configurations, with small 7-point mosaics 
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in the compact D and E configurations; the resulting 
images are sensitive to spatial scales from 10 2 to 10 4 AU. 
Our high resolution maps are also sensitive to binaries 
with separations > 250 AU. In only one source is there 
clear evidence for multiple components within 2000 AU, 
suggesting a tentative multiplicity fraction of only 10% 
for very deeply embedded protostars. 

For un-blended sources, our CARMA observations re- 
cover 50 — 95% of the single dish Bolocam flux, assuming 
a dust spectral index of /3 = 1.5. In all nine observed 
sources, 230 GHz emission is seen on both large scales 
(a few 1000 AU), indicating an extended envelope, and 
small scales (down to 170 AU), indicating a compact disk 
component. Three of the nine Serpens Class sources 
were not observed here due to non-detections or marginal 
detections in preliminary 110 GHz CARMA data, and 
these sources likely do not have a significant disk. Thus, 
we see evidence for circum-protostellar disks in 6/9 Class 
sources in Serpens. 

We make a first estimate of the disk mass using the 
visibility amplitude at 50 kA, assuming an optically thin, 
unresolved disk with Tp = 30 K and a distance to Ser- 
pens of 415 pc (|Dzib et al.ll2010f ). Resulting disk mass 
estimates range from 0.04 M to 1.7 M , with a mean 
disk mass of 0.2 M when the three non-detections are 
included. Errors in disk masses may be ±50%, depend- 
ing on the envelope contribution at 50 kA, the size of the 
disk, and the presence of optically thick emission. For 
a distance to Serpens of 260 pc, as assumed in our pre- 



vious work, masses would deccrease by a factor of 2.5. 
Our mean Class disk mass is similar to that of the 
la rger, but less wel l defin ed, Class and Class I sample 
of iJdrgensen et all (|2009f ). 

Our high disk detection rate in a complete sample 
suggests that circum-protostellar disks are quite com- 
mon even in the youngest protostars, with approximately 
2/3 of embedded protostars in Serpens having detectable 
(M disk > 0.02 M ) disks by 0.2 Myr. 

Future work comparing the high-resolution millime- 
ter data, as well as spectral information, with radiative 
transfer models will provide a more robust determination 
of the disk masses. 
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APPENDIX 
PRELIMINARY 110 GHZ DATA 

Preliminary v = 110 GHz (A = 2.7 mm) observations of the Bolocam 1.1 mm cores Ser-Bolo 8 (associated with 
mid-infrared identified protostars Ser-emb 3, emb 7, and emb 9) and Ser-Bolo 24 (associated with Ser-emb 2 and 
emb 15) were completed with CARMA. The 6.1-m and 10.4-m antennas were used to obtain 110 GHz continuum 
observations in the D and E antenna configurations between 2007 April 20 and 2007 September 04. 

Three correlator bands were configured for continuum observations with 468.75 MHz bandwidth, for a total band- 
width of 2.82 GHz. 1751+096 was observed approximately every 20 minutes for complex gain calibration. Absolute 
flux calibration was accomplished using 5 minute observations of MWC 349. The overall calibration uncertainty is 
approximately ±20%, from the reproducibility of the phase calibrator flux on nearby day s. The passband ca librator, 
1751+096, was observed for 15 minutes during each set of observations, and either optical (jCorder et al. 2010) or radio 
pointing was performed every one to three hours. Data reduction was completed as described for the 230 GHz data 
in Section [U 

The resulting 110 GHz maps are shown in Figure [TUJ The rms noise in both maps is ~ 1.5 mJy, and the synthesized 
beam sizes are 5" — 6" '. Ser-emb 7 and emb 15 are both strongly detected. Ser-emb 3 and emb 9 are undetected at 
the 3cr = 4.5 mJy level, although there is some 2a emission at the position of both sources. Ser-emb 2 is marginally 
detected at the 3<r level. 

Any detected flux in the compact D and E configurations may well include envelope emission, so we adopt an upper 
limit of 4.5 mJy for the disk flux of all three sources. We convert this flux upper limit to a disk mass upper limit using 
equation [TJ where K^mm = 0-003 c m 2 g - 1 is the dust opacity at 110 GHz. The opacity is extrapolated from Table 1 
column 6 of Osse nkopf fc Henningl (jl994j ) assuming n oc v& and /? = 1.5. 

This yields a disk mass upper limit of 0.02 M for Ser-emb 2, emb 3, and emb 9. While the preliminary 110 GHz 
maps are less sensitive than the subsequent 230 GHz observations, non-detections place a fairly stringent upper limit 
on the total envelope mass, and thus on the disk mass. 
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